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Abstract 

A generalized inverse scattering method has been apphed to the 
hnear problem associated with the coupled higher order nonlinear 
schrodinger equation to obtain it's iV-soliton solution. An infinite 
number of conserved quantities have been obtained by solving a set of 
coupled Riccati equations. It has been shown that the coupled system 
admits two different class of solutions, characterised by the number of 
local maxima of amplitude of the soliton. 
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1 Introduction 

Optical soliton since it's discovery occupies a distinguished position in the 

nonlinear optics research and it is regarded as the next generation carrier in 
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an all-optical communication system. Because of their remarkable stability 
solitons are capable of propagating long distance without attenuation [1, 2]. 
Theoretical study of optical sohton began as early as in seventies with the 
landmark discovery by Hasegawa and Tappert [3] . They had shown that the 
nonlinear Schrodinger equation(NLSE) studied by Zakharov and Shabat [4] 
(but with space and time interchanged) is the appropriate equation for de- 
scribing the propagation of pico-second optical pulses in optical fibers. Later 
on MoUenauer etal. [5] successfully demonstrated the transmission of pico- 
second optical solitons through a monomode fiber. Recently Kodama and 
Hasegawa [6] proposed a modified NLS model, known as the higher order 
nonlinear Schrodinger equation (HNLSE), which is suitable for the propa- 
gation of femto-second optical pulse. Untill now only two integrable HNLS 
equations [7, 8] which has soliton solutions are known. The dynamics of 
HNLSE not only takes care of dispersion loss but also takes care of the prop- 
agation loss as the optical sohton propagate along the fiber. This is due to 
the fact that the stimulated Raman scattering effect, which compensates the 
propagation loss, already exists in the spectrum of the HNLS equation. As 
a consequence a short pulse can propagate as a soliton for a long distance 
without distortion. Therefore the importance of the study of HNLS system 
and to find their soliton solution is unquestionable. The soliton solution may 
be obtained by several methods hke Hirota's bilinear method [27, 28, 29], Lie 
group theory [30, 31, 32] and inverse scattering method (1ST) [4, 33, 34]. The 
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1ST is the most elegant and one of the most widely used technique, which 
eventually proves the complete integrability of the system. 

It has been known that a good number of soliton equations in nonlinear optics 
have multifield generalizations. The coupled nonlinear Schrodinger equation 
proposed by Manakov [9], is one such example. It describes a system with 
two interacting optical fields with different states of polarizations. It has an 
important application in communication systems using optical solitons. Such 
systems have many other important applications in all optical computations 
[10, 11, 12, 13], in photorefractive crystals [14, 15]. Recently the existence of 
multicomponent solitons has been discovered experimentally also [16] . Inter- 
estingly the HNLS equation also has it's multifield generahzation and they 
have been studied in different physical contexts [17, 18, 19, 20]. The cou- 
pled HNLS (CHNLS) equations are particularly important for describing the 
dynamics of ultrashort optical pulse in a system involving two or more in- 
teracting optical fields. Recently a CHNLS equation, which incorporates the 
effect of third order dispersion, kerr nonlinearity and stimulated Raman scat- 
tering was proposed by Nakkeeran et. al., [21, 22]. The equation considered 
in [21, 22] is given by 

• n n n 

^ j=l 3=1 j=l 

(1) 

for A; = 1, 2, • • - n. It describes the evolution of a complex vector field E 
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with n components, where n is a finite integer. The suffices z and r denote 
the normahsed space and time derivatives respectively. Ek and E'^ respec- 
tively represent the amplitude of the kth component of a slowly varying field 
and it's complex conjugate, e is an independent parameter. They have also 
shown the integrability of (1) through the existence of Lax pair and obtained 
the one soliton solution (ISS) of (1) through Backlund transformation. Ex- 
istence of Lax pair although indicates the integrability of the system but a 
conclusive sense of integrability is achieved if the system possess an infinite 
number of conserved charges [26]. It is important to note that the equation 
(1) is the multifield generalization of the Sasa-Satsuma equation [7]. Inter- 
estingly for Sasa-Satsuma (scalar) equation two different ISS were reported 
in [7, 23, 24, 25]. One is the convensional sech solution having a single peak 
and the other is a complex combination of sech function having two peaks. 
In view of the connection between the Sasa-Satsuma equation and it's multi- 
field generalization it is expected that the latter equation should also admit 
both type of soliton solutions. It should be mentioned that the ISS with 
sech envelope function has been obtained in [21, 22]. However, it remained 
to show that (1) admit a more general class of solutions. The existence of the 
higher order solitons, which are important to study the soliton interactions 
also remained unexplored. A possible reason for this may be the difficulty 
involved in solving the linear problem associated with the system. It is im- 
portant to note that for a completly integrable system it is necessary that 
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the system possess A^-soliton solution. 



Our objective in this paper is to obtain the whole hierarchy of the conserved 
charges in a systemetic way to establish integrability of (1) conclusively. We 
complete the investigation of higher order soliton solutions for (1) by deriving 
the exact A^-sohton solutions for the system. We invoke the inverse scattering 
method for an (2n+l x2n+l) dimensional eigenvalue problem and obtain the 
A^-soliton solution by solving the set of (2n + 1) Gelfand-Levitan-Marchenko 
(GLM) equations. Subsequently we discuss some important characteristics 
of the ISS. 

The paper is organized in the following sequence. In section 2 we use a set of 
variable transformations to cast the CHNLS equation into a form, suitable 
for the inverse scattering transform. Subsequently we study the linear prob- 
lem associated with the modified system. In section 3 we show the existence 
of infinite number of conserved quantities by solving a set of coupled Riccati 
equations. The 1ST scheme and the properties of the scattering data are 
studied and the generahzed GLM equations are derived in the section 4. The 
exact A^-soliton solution is obtained in section 5. In Section 6 two class of so- 
lutions of the CHNLS equation are discussed. Section 7 is the concluding one. 
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2 Eigenvalue problem 

In order to investigate (1) we use the following change of variables and a 
Galelian transformation, 

qk{x,t) = Ekiz,T)e-^e^^-^e\ (2) 
t = z, (3) 

X = T — (4) 

Thus we get a set of complex modified K-dV equations (CMKDV) 

n n 

We use X and f in (5) and in the subsequent expressions to denote the deriva- 
tives with respect to x and t respectively. The associated spectral problem 
of (5) may be represented by a pair of linear eigenvalue equations. 

a,* = L(a;,t,A)* (6) 

= M(a;,t,A)^ (7) 

where ^ is the {2n + 1) dimensional auxiliary field and A is the time inde- 
pendent spectral parameter. The pair of matrices (L, M) is called the Lax 
pair. The explicit form of the Lax pair has been given in [21, 22]. With slight 
modifications (see eq. 11 in [22]) we express the Lax pair (L, M ) by using a 
pair of matrices (S,/!), where S is a c-no. diagonal matrix and the matrix 
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A{x, t) is a potential function of the eigenvalue problem and consists of the 
dynamical fields, qk{x,t) and ql{x,t) only, 

2n 

S = X/ ^'^f' ~ ^2n+l 2n+l 
k=l 

n 

A{x,t) = ^(qk{x,t)e2k-l 2n+l 
k=l 

+ ql{x, t)e2k 2n+i - ql{x, t)e2n+i 2k-i - qk(x, t)e2n+l 2k) (8) 

where, e^j is an {2n + 1 x 2n + 1) dimensional matrix whose only {kj)th 
element is unity, the rest elements being zero. By using the properties of E 
and A 

viz, 

EA + AE^O, = 1 

we write the Lax matrices in a simplified form, 
L = -iAE + A 

M = e{-A^^ + 2iXEA^ - {AA^ - A^A) - 2i\Y.A^ + 2A^ + AX^A - 4iA^E) 

(9) 

The Lax pairs given by (9) is also valid upto an additive constant, since a 
constant commutes with all other matrices. We assume the constant to be 
zero for simplicity. The consistency of the Lax equations (6, 7) gives the 
nonlinear evolution equation for the matrix A, 

At + e{A,^^-?,{A,A^ + ?,A^A,)) = Q (10) 
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which is nothing but the matrix form of the equation (5). Although the ex- 
istence of a Lax pair is itself a sign of the integrability of (5) vis. a. vis. (1), 
it remains to show that the system also possess infinite number of conserved 
quantities and admit A^-soliton solutions. 



3 Conserved quantities 

In order to obtain the conserved quantities first we derive the associated Ric- 
cati equations. To this aim we write the Lax equation (6) in the component 
form. For the first 2n components of ^, the equation (6) can be written in 
the form 

*2fe-l X = -iA*2fc-l + ?fe*2n+l 

*2fex = -iX^2k-l + ql^2n+l (11) 



for A; = 1, 2, • • • , n. But the {2n + l)-th component has a different form 

n 

*2n+l X = iX'^2n+l " ^(?fe*2fe-l + gfe*2fe) (12) 

k=l 

Following now a similar procedure as in [24] we write, 

r. = (13) 
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for k — 1, 2, • • ■2n. By using (11, 12, 13) we may obtain a set of first order 
differential equations, 

n 

r2jk-i + 2iXr2k-i - ^{qlT2k-i^2j-i + qk^2k-i^2j) - Qk = (14) 

n 

T2k + 2zxr2k-Y.i^kr2kr2j-i + qir2kr2j)-ql = o (is) 



wfiicfi are known as tfie Riccati equations. Tfie solution of the equations are 
related to the conserved quantities {a2n+i 2n+i) in the following way 



/oo " 
dxJ2{qkr2k-i + qtr2k) (16) 
-■^ fc=i 



lna2n+i 2n+i(A) = - iXx 

We will see in the section 4 that Q;2n+i 2n+i is the diagonal element of the 
scattering matrix and it is time independent. 

In order to solve Riccati equations (14) and (15) we assume r2jk_i and 
r2k in the form 

oo 

^2k-i = Y.Cf-\x)W> (17) 

p=0 
oo 

^2k = T.Cl\x)\^ (18) 

p=0 



Substituting (17,18 ) in (14) we get 

Cf-^ = 0, Cf-^^^. (19) 

n p+1 
1=1 m=0 
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similarly substituting (17,18) in (15) we get 

= 0, Cl' = % (21) 

n p+l 

2iClX, = -(C|J J. + E E + (22) 

1=1 m=0 



Now the infinite number of Hamiltonians (conserved quantities) may explic- 
itly be determined in terms of the q^i ql and their derivatives by expanding 
a2n+i 2n+i in the form 

ln{a2n+l 2n+l) = 71^-—^//^ ^ (23) 

The first few conserved quantities are given by 

Hi ^ - [ dxJ^QUk (24) 
fe=i 

1 . n n 

H, = _ /dx(2 5:|gfc|Y-^g,%5fe,) (25) 
^ •' k=i k=i 

-1 „ n n n 

H, = -/dxE9^..g...+2(Ek^r)'-(E(k^rw' 

" •' k=l k=l k=l 

n n n n 

- Y.l^kl'^Y.llxQkx-Y.^kQkxY.lixQix (26) 

k=l k=l k=l 1=1 



Notice that the hamiltonians with odd indices only survive. The even indexed 
hamiltonians become trivial. It is easy to see using the equation of motion 
that Hi, H3, are indeed the constants of motion. If we choose the field q, a 
scalar field then we are able to show that Hi,H^, H5 reduce to the conserved 
quantities for the scalar Sasa-Satsuma equation [24] . 

11 



4 Gelfand-Levitan-Marchenko equations 

We now generalise the inverse scattering method for the (2n + 1 x 2n + 1) 
dimensional Lax operators (9). The generalization however, is a nontrivial 
one and crucially depends on the scattering data matrix, we have broadly 
followed the treatment of Manakov, developed in the context of 3 x 3 Lax 
operators [9]. 

In order to formulate the scattering problem we assume that the family of 
Jost functions $(fc=i>2, -2n+i) ^^^^ ^(fe=i,2,...2n+i) ga^igfy the following 

boundary conditions for real values of A, 



e^e-*"^ (27) 



for A; = 1,2... 2n, but the (2n+l)-th component satisfies a different boundary 
condition 



$(2n+l) 



e2n+ie'^ (28) 



Similarly, other set of Jost functions satisfy the boundary conditions. 



e^e-^^^ (29) 



for A; = 1, 2 ... 2n and the (2n + l)-th component satisfies. 



^(2n+l) 



e2n+ie'^^ (30) 



In the equations (27-30) e^'s are the basis vectors in an {2n + l)-dimensional 
vector space. Note that the set of jost functions (27-30) also satisfy the 
orthogonality condition. That is, 
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for k,j = 1, 2, . . . 2n + 1. Since vectors ^''^'^^ form a complete set of solutions 
of (6) hence, 

2n+l 

$(^)(a;,A)= J2'^kj{X)^^'\x,X) (32) 

where Q;fej(A) is the {kjY^ element of scattering data matrix (det[Q;fcj]=l). 
Using (31) and (32), akj is expressed in the form 

ak,{X) = ¥^^\x,X)^^''\x,X) (33) 

It is interesting to see that using the unitary property of [akj] we can write 
a^j as the cofactor of the elements of the matrix [a^j], that is, 

<+i k = (-l)'"+'+'rfet[a2n+i k] (34) 

where [(52n+i k] is an (2n x 2n) dimensional matrix, constructed from the 
matrix [ajkj]2n+i 2n+i with (2n + l)-th row and k-th column being omitted. 
Now by using (32) and (34) we obtain the following useful relations among 
the jost functions, 

1 ^"^ rv* (\) 

'^2n+l 2n+l['^) j=l '^2n+l 2n+l['^) 

(35) 

for A; = 1, 2, • • • 2n. The (2n + l)-th jost function satisfy another relation, 

1 1 

^(2n+i)g-ax ^ ^(2n+i)^-iXx ^ ^ ^^^^^ ^-^^e-^^^ (36) 



Oi2n+l 2n+l Oi2n+l 2n+l ,=i 
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Notice that in deriving (35, 36) we have used the following property of scat- 
tering matrix, 

2n 

^2n+l k^ij = J2[^2n+l k\il(Adj[a2n+l k\)lj (37) 
1=1 

In order to obtain the complete analytic behaviour of the jost functions vis 
a vis scattering data the domain of A is extended to complex plane. It can 
be shown that the functions $('^)e*'^^ for A; = 1, 2, • • • 2n and ^(2n-+i)g-«Aa; 
analytically continued into the upper half- plane {ImX > 0) whereas ^*('=)e*-^^ 
for k — 1,2, ■ ■ ■2n and $*{2n+i)g-«Aa; g^^g analytically continued into the lower 
half-plane. Consequently the scattering element, Q!2n+i 2n+iW 
ments of the matrix [Q;2n+i 2n+ii^)] are analytic in the upper half-plane and 
0!2n+i 2n+ii^) ^-nd all elements of the matrix [a2n+i 2n+iW] ^-re analj^ic in 
the lower half-plane. It is important to note that the bound states of the 
eigenvalue equation (6) correspond to zeros of a2n+i 2n+iW in the lower half- 
plane. We assume that the bound states are located at {ImX > ) for 
j — 1,2, ■ ■ ■ N, where the jost function becomes, 

2n 

$(^"+i)(a;. A,*) = Y: C^i+i rJ^^'^Kx. (38) 

m=l 

In (38), C2n+i m represent the value of the scattering parameter a2n+i m at 
the position of the pole. 

The time dependence of the scattering data may be easily obtained from the 
asymptotic limit of (7), which gives the following time dependence of the 
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scattering data. 

a2n+i kit) = a2n+i fe(0)e-''^^'* (39) 

a2n+l 2n+l(^) = "2n+l 2n+l(0), (40) 

^iii kit) = .(0)e-«-^^^* (41) 

C2n+1 2n+l(^) = C's^+i 2n+l(0) (^2) 

In order to derive the GLM equation we consider an integral representation 
of the Jost functions 

*(^)(x,A) = 6^6-*^^+ / dyK^^\x,y)e-'^y (43) 

Jx 

with j = 1,2, . . . 2n, while the (2n + 1)*^ Jost function is considered as 

roc 

(a;, A) = e2n+ie*^^ + / d^/K^^^+i) {x, y)e'^y . (44) 

where, the kernels K^^^ and K^^'*"'"^) are (2n+ 1) dimensional column vectors, 
which may be written explicitly in the component form as 

2n 

^^'\r.y)= K^ri\r.y)em (45) 

m=l 

2n+l 

^^'-^'\x,y)= E^i"^'n^,y)en. (46) 

m=l 

Multiplying (36) with ^ /_!^ e-'^'HX, [y > x) and using (43,44,38) togather 
with the analytic properties of the associated scattering data we obtain the 
desired GLM equation for the kernet 



2n 2n 



^"^ POO 

K(^"+^)(a;, y) + Y. ^P^pi^ + y) + 12 dsk^^\x, s)F,is + y) = (47) 
p=i p=i ''^ 
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where, 

The ' over a2n+i 2n+i denotes derivative with respect to A. 

The other integral equations for kernel is obtained from (36) in a way 
similar to that (47) and the resultant equations are 

K^\x, y) + e2n+iF*{x + y) + T dsk^^^-^'\x, s)F*{s + y) = (49) 

Jx 

for p = 1, 2, • • • 2n In deriving (49) we have used the identity, 

2n 

Qn+l m = «2n+l mi^j) = IZ[«2n+l mi>^j)]kiiAdj[a2n+l mi>^j)])il (50) 

i=l 

The set of equations (47,49) may be called generalized GLM equations. Sub- 
stituting (46,49) in (47) we get the GLM equation for the p-th component of 
of the kernel A;^^""'"^-*, which is given by, 

^ , roc 

{x, z)+Fp{x + z)+J2( dsKf "+^^ (x, s) 

m=l ■'■■^ 

POO poo V 

+ / dyF2m-i{y + z)F;^_^{y + s) / dyF^Jy + z)F2m{y + s)) = 

(51) 

5 N-soliton Solution 

To obtain a closed form solution of the GLM equation (51) we assume 
{a2n+\ p(A) = 0). This is justified because our primary interest is to obtain 
the soliton solutions, which is obtained for the reflectionless potential. We 
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also note that, substituting (44) in (6)we get a relation between the Kernels 
of the integral equations (51) and the 'potential' of the eigenvalue equation 

(6), 

q,,.,{x) = -2Kf-f\x,x) (52) 

qUx)^-2Kfr'\x.x) (53) 

for i = 1, 2, • • • n. To obtain the general A'"-soliton solution diagonal element 
of the scattering matrix, Q;2n+i 2n+i (A) is cosidered to have A'"-pairs of zeros 
located symmetrically about the imaginary axes in the lower-half plane. That 
is. 

Note that unlike the CNLSE, where a zero corresponds to a single soliton, 
in the CSSE a pair of zeros corresponds to a single sohton. Finally we solve 
the set of GLM equations, by assuming that kernel of the integral equations 
are of the form, 

TV 

(x, z) = Y, Rpj{x, t)e-^> + Spj{x, t)e^^-^ (55) 
and we obtain the /^-sohton solution of (5) for the A;*'* component of the field, 

2N 

Qkiz, t) = -2 ^(BC-^),,e-^^*^ (56) 

B and C in (56) are respectively 2n x 2N and 2A'" x 2A'" matrices whose 
elements b^j, c^j are given by. 



Pfe.e-^"^ 1 < J < iV 
N +1 < j <2N ^ 
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(57) 



and 

2n 



(Ar.A^,.)(A..A;,.) ^ l<m<iV, 

^ / , Pkm Pk] 6 -' 

T. '^' ixt.>^iX.^>^ V N + l<m<2N, 

j=i 



where, 



with the constraints 



P>,=. /"^'''f (59) 

Q; 2n+l 2n+l{'^j) 



A,+iv = -A* (60) 

Pkj+N = Pfej (61) 

for J = 1, ■ ■ ■ , -?V. The A?"- sohton solution for the equation (1) is obtained 
from (57) by using the inverse variable transformations (2,3,4). 

6 One soliton solution 

The one soliton solution for (5) is obtained by choosing A?" = 1 in (56). It 
implies from (54), that a2n+i 2n+i(Aj) has a pair of zeros located symmetri- 
cally about the imaginary line in the lower half-plane. We assign them as 
(— Ai, A^), where Ai — {—^ + ir])/2 with ^, 77 > 0. By reinstating the transfor- 
mations (2,3,4) in (56) we obtain the k-th component of one soliton solution 
for (1), which is given as. 
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where 



Ek{z,T) = 



27] pki e"- {e + c e ^) 



In 



X]pj-i(|c|-ie2^ + |c|e-2^ + 2|c|) 



(62) 



12e2' 
1 



1 



12e2 108e3^ 



c = 



1 

1-,- 

TjC* 



(63) 
(64) 
(65) 

(66) 



Notice that each component of the sohton solution (62) is defined completely 

2n 

by a set of four parameters viz, r],^, pki, and ^P/i- It represents an envelope 
wave moving with a group velocity, e(?7^ — 3^^ + j^) undergoing internal 
oscillation. Interestingly the group velocity depends on both real part (77) 
and imaginary part (^) of Ai. If we specialize to the scalar limit of the 
CHNLS equation the solution (62) reduces to the one soliton solution of the 
Sasa-Satsuma equation (see Eqs 42. in [7]). In order to investigate the shape 
of the pulse we take the derivative of |-Efcp with respect to r. This gives the 
following conditions for maxima of l^'fel^, 

e'^ = |c|2 - 2 ± ^(|cP - 2)2 - |cP = (67) 
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It is clear from (67) that for |c|, such that 1 < |c| < 2 the solution has a 
single peak. Interestingly however, for |c| > 2, there are two values (real) 
for e^"^, that is (62) has two maxima. This corresponds to two peaks. As 
|c| increases further the two peaks gradually shifts apart from each other. 
Finally at |c| — > oo, that is, when the two zeros of a2n+i 2n+i merge on the 
imaginary line in the lower half-plane the solution (62) then reduces to. 



sech{r]T - er]^z + —z)e'^^- (68) 



2n ^' ' 12e 



It is important to note that the phase factor arises in (68) is purely from 
the variable transformations (2,3, 4). The solution (68) represents an wave 
moving with a group velocity, {erf + -^). Unlike the earlier case (62), the 
group velocity depends only on the real part 77 of Ai. This is the solution 
reported in [21, 22]. It is important to note that the shape of the solitons 
remains invariant with respect to space and time for all values of |c|. The 
two class of solutions obtained for ISS may be extended straightforwardly 
for N- soliton solutions. 



7 Conclusion 

In this paper we have studied the CHNLS equation by applying a gener- 
alized inverse scattering method developed to solve the (2n -|- 1 x 2n -|- 1) 
dimensional linear problem associated with (5) vis a vis (1). We have shown 
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the integrability of the system by showing the existence of infinite number 
of conserved quantities. The A'"-sohton solutions for the system have been 
obtained by solving a set of generalized GLM equation. We have shown two 
different class of solutions by considering the zero's of the diagonal element 
of the scattering data on the imaginary line and a pair of zero's lying sym- 
metrically about the imaginary line in the lower-half plane. By a suitably 
defined parameter we have shown how the double-peak soliton reduces to the 
single peak sohton. The results, we have obtaind predicts that CHNLS equa- 
tion allows dispersionless propagation of the ultrashort optical soliton in the 
shape of single hump pulse or double hump pulse. This may have interesting 
consequences in the propagation of optical solitons through nonlinear fiber. 
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